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Two systems are home-metric if they are indistinguishable by diffraction. We first make a dis- 
tinction between Bragg and diffuse scattering homometry, and show that in the last case, coherent 
diffraction can allow the diffraction diagrams to be differentiated. The study of the Rudin-Shapiro 
sequence, homometric to random sequences, allows one to manipulate independently two-point and 
four-point correlation functions, and to show their effect on the statistics of speckle patterns. Con- 
sequences for the study of real materials are discussed. 
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INTRODUCTION 

The possibility to shape coherent X-ray beams from 
synchrotron light sources [1] and to get naturally coher- 
ent beams from X-ray free electron lasers (XFEL) Q, 
has revolutionized the way X-ray diffraction experiments 
are performed and analyzed. One of the most fascinat- 
ing property of coherent diffraction is the possibility to 
measure speckle patterns [3], which are much more in- 
formative than the diffuse scattering obtained by classi- 
cal diffraction. Together with the development of novel 
sources, phase retrieval algorithms have also emerged, al- 
lowing the reconstruction of the diffracting objects under 
certain conditions 0, Hj]. However, the reconstruction of 
a structure is not always possible nor necessary to study 
the physics of materials. For example, measuring corre- 
lation lengths close to phase transitions [6[ or slow dy- 
namics with X-ray photon correlation spectroscopy [H, uj 
does not require the full reconstruction of the system un- 
der study. 

The purpose of this letter is to show that statistical 
analysis of speckle patterns can yield information on or- 
ders hidden to conventional X-ray analysis. In this re- 
spect, we are in line with recent works showing that four- 
point intensity cross-correlation of speckle patterns can 
uncover " hidden symmetries" present in colloidal glasses 
[1, Q or magnetic systems [loj . 

Our approach uses the concept of homometry, i.e. the 
property of different systems to exhibit same diffraction 
patterns. Separating out the scattered intensity expres- 
sion into three terms allows one to show that homometry 
can occur at different levels. We then put the emphasis 
on diffuse scattering homometry, that we discuss with the 
help of the well-known Rudin-Shapiro sequence 11 , HI- 



COHERENT DIFFRACTION 

Let us first give a general expression of the intensity 
scattered at scattering vector q, by a ID periodic TV- 
site lattice decorated by two atoms A and £?, of scat- 
tering factor Ja and in proportion x and 1 — x re- 



spectively. Generalization to 2D, 3D, multi- atomic basis, 
displacement disorder, or disorder of the second kind [HI 
is straighforward. Following ref. [13], the diffracted in- 
tensity is given by the formulae: 

I(q) = E /»/»'e i9( "'- n) = E E fnfn +m e i<im . (1) 



The ensemble average of the cross-product (fofm) is then 
introduced: 
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where N m is the m-dependant number of terms of the 
sum J^ n . The A m term, usually negleted in textbooks, 
is due to finite-size fluctuations of the spatial average 
with respect to the ensemble one. 

Further introduction of Af m = f m — (f) allows one to 
get the three components of kinematic diffraction: 
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The first term gives the intensity of the Bragg reflec- 
tions, and the fringes due to finite size effects. For a 
crystal of N cells of structure factor F(q), it can be writ- 
ten as: 



I B (q) = \(F(q))\ 



2 sin 2 q(N + l)/2 
sin 2 q/2 



(4) 



In practice, the fringes given by the sine functions are 
only visible with a coherent beam illumination (for con- 
ditions of observation and examples see |3|)- 

The second term is the diffuse scattering intensity, 
which only depends on pair correlation function (CF). 
For random disorder, it reduces to the well-known Laue 
formula: 



I(q)=Nx(l-x)(f A -f B ) 2 , 



(5) 
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FIG. 1. (left) (40 x 40) lattices of +1 and -1 in equal pro- 
portion and (right) diffraction patterns in the first Brillouin 
zone (log scale) for a) a triplet SRO lattice [3] b) a random 
lattice. Reciprocal nodes are in the corners. 
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where random refers to the vanishing of the pair CF (i.e. 
(A/ A/ m ) = (A/ )(A/ m ) = for m ? 0.) 

The third term gives rise to speckles. Like fringes, 
speckles only exist if the incident beam is coherent 
enough, and if the system does not explore too many 
configurations during acquisition time T (non-ergodicity 
condition (A m )x ^ 0). Interestingly enough, working 
out of coherent conditions has the effect of averaging out 
A m , which yields ensemble averaged quantities. 



HOMOMETRY 

Homometry - etymologically same distance - is a word 
coined by A. Patterson [15|, [l6|, to describe the prop- 
erty of non-congruent sets of points to possess the same 
pair distances (or the same difference sets) [I?}- Homo- 
metric sets have thus the same diffraction pattern, as 
Eq. (pQ) demonstrates. A simple example of homome- 
try is given by the two sets S = {0, 1,4, 10, 12, 17} and 
S' = {0, 1, 8, 11, 13, 17} Indeed, the structure fac- 

tors F(q) of the two sets have the same magnitude for 
all q-vectors, but not the same phase. Hence, the lost 
of the phase makes these sets indistinguisable by X-ray 
diffraction. 

However, because solid state physics deals with ma- 
terials, the above definition turns out to be too restric- 
tive. Eqs (j3j) allows to distinguish between Bragg (B) ho- 
mometry, diffuse scattering (D) homometry and coherent 
diffraction (C) homometry. 



BRAGG HOMOMETRY 



FIG. 2. Diffraction patterns and associated intensity varia- 
tions from ((iV, M) = (64,512)) a) RS and b) BS sequences. 
The vertical broadening of the patterns are for visual conve- 
nience. 



sented in [16j, of unit cells size equal to 8 and atomic 
positions given by H = {0,3,4,5} and H' = {0,4,5,7}. 
Structure factors, readily calculated as: 



Fh(q) = 1 + 2cosg + cos4g 
Fw(q) = 2(cos2<? + cos4<?), 



(6) 
(7) 



B-homometry decribes crystals with different basis but 
same Bragg intensities [15l Il6j. To illustrate that, let us 
consider the examples of ID homo metric crystals pre- 



have the same amplitude squared at the Bragg positions 
q = 27r/i/8 but not out of Bragg positions. Eq ((4]) shows 
that the fringes intensity, revealed by coherent diffrac- 
tion, gives out-of-Bragg values of \F(q)\ 2 which, at least 
in theory, allows to distinguish H and H' and solves the 
B-homometry issue. This is well-known and corresponds 
to the oversampling requirement of the phase retrieval 
algorithms [18|, LL9| . It is clear however that if the atomic 
basis is homo metric itself, like e.g. in crystals with S or 
S' basis, the problem cannot be solved, coherence or not. 



DIFFUSE SCATTERING HOMOMETRY 

Surprising examples of D-homometry were designed 
by Welberry et al. [l3, Eoj. They consist in substitu- 
tionnally disordered lattices with triplet (or quadruplet) 
short range ordered (SRO) CF [21(, but zero two-point 
correlations (Fig. [2ji). Diffraction diagram of these lat- 
tices present the same Bragg and diffuse scattering inten- 
sity [14[. But, as shown in Fig. [2j their speckle patterns 
are different, which breaks the D-homometry. 

A more tractable (though subtle) example of D- 
homometry is provided by the geometrically ordered 
(GO) [HI Rudin-Shapiro (RS) sequence Q, whose 
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FIG. 3. Probability densities P(I) of diffraction patterns of 
(N,M) = (4096,4096 x 128): a) BS (green) and RS sequence 
(blue), b) sequences obtained by Bernoullization c) #4-SRO 
sequences for f - N/37 (p = 0.99), N/13 (p = 0.95) and 
N/3 (p = 0.999). The line indicates the random negative 
exponential law. 



generic term a n can be written [12[: 

__ j a n for I = 0, 1 

a4n+z ~ 1 {-l) n+l a n for Z = 2,3 



with ctq = 1. (8) 



This sequence has become famous [ll|, l22M24j because, 
thought GO, it is D-homometric to randomly distributed 
sequences (sometimes called Bernoulli sequences (BS) 
j24j), with the same ANx(l — x) diffuse scattering in- 
tensity (Eq. [5]). In other words, its two-point CF 
92(n) — o-o&n is zero for n ^ (spatial average). 

Let us now consider the coherent diffraction of RS and 
BS sequences. In what follows, we present Fast Fourier 
Transforms (FFT) computations of sequences of length 
iV, zero-padded up to a value M ^> N to clearly see the 
speckles. For the RS sequence, because x / 0.5 and de- 
pends on its length iV, we have always susbstr acted the 
average value 2x — 1 to all the terms in order to get rid 
of the Bragg intensites. BS of 1 and -1 in equal pro- 
portion were computed with the python pseudo-random 
number generator. The squared value of the FFT I(q) 
are normalized by 4Nx(l — x) in order to get I(q) = 1. 

Figure [2] shows the diffraction patterns of RS and BS 
sequences in the first Brillouin zone (0 < q < 2tt). Inspec- 
tion of these patterns shows that, although their average 
value is the same, the speckles repartition is remarkably 
different. In particular, it is clear that the BS pattern 
exhibits much more spikes, while the RS pattern is more 
homogenous and regular (no low-intensity speckles more 




FIG. 4. a) Computed g^n) using periodic boundary condi- 
tions for N = 256 RS sequence (blue) and a BS sequence 
(green), b) Magnitude of RS g^q) as a function of h — q/2n. 



regularly spaced). To quantify this observation, we stud- 
ied the statistics of both speckle patterns by calculating 
their probability density of intensity P(I). It is known 
that for a random media (see e.g. [25]), the intensity dis- 
tribution has a negative exponential distribution given 
by: 



P(I) = jexp-j, 



(9) 



which in our case reduces to P(I) = exp(— I). 

Figure [3^) shows the probability densities P(I)rs and 
P(I)bs for the RS and BS sequences. While P(I) B s 
follows quite well the negative exponential law, as ex- 
pected, it is not the case for P(I)rs- Though the pre- 
cision of P(I)rs depends on N, it is well approximated 
by the step function P(I < 2) = 0.5. This statistics, 
which means that intensities lower than 2 occur with the 
same probability, explains the homogeneous aspect of the 
diffraction pattern. The reason for this unusual statistics 
is not clear, but undoubtedly comes from the GO na- 
ture of the RS sequence. The presence of order, invisible 
through diffuse scattering, is revealed by the statistics 
of the speckle pattern, breaking the D-homometry in a 
quantifiable way. 

In order to test the robustness of the Prs (I) behavior 
with respect to disorder, we first quantify the degree of 
order of the RS by one of its quadruplet CF: 



g±(n) = ctq g\ (j n <r n _|_ i . 



(10) 



Indeed, we found numerically that at variance with the 
BS, g 4 (n) is LRO for the RS sequence (Fig. Ht)^- This 
is confirmed by the behavior of its FFT g±(q) (Figure HJd), 
which exhibits well defined peaks indexed by the basis 
vectors {j^\hi G j-1, G N} [l?} characteristic of 
limit-periodic functions [28| . By analogy with two-point 
orders, we define 774 = g^ir /2) /N 2 = 1/4 as the order 
parameter of this sequence. 
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We first decreased 774 while keeping g 4 -LRO and #2- 
disorder by using the " Bernoullization" procedure as de- 
fined in [12J. It consists in changing the sign of each <j n 
with probability p, in order to build sequences interme- 
diate between the pure RS (p = 0,1) and BS (p = 0.5) 
sequences. The order parameter 774 was found to vary as 
m (p) ~ (1 - 2p) 4 /4. 

Typical density probabilities density shown in Fig. 
[3b), exhibit a continuous evolution as a function of p. 
The step- function behavior is rapidly lost as p — > 0.5, 
with the best sensitivity close to the small intensity val- 
ues P(I = 0) = Po- Simulations show that Po(p) closely 
follows 1 — 2774 (p), and that sizeable deviation from the 
normal law starts from p > 0.8. 

Effect of #4-SRO on the speckle pattern was studied 
by shifting the sequence by two lattice periods at certain 
points, randomly selected with probability p. This en- 
sures to achieve #4- SRO, clearly seen by the broadening 
of the #4(71-/2) peak, while keeping the g2{n) correlation 
to zero. Figure [3fc) shows the results for different p val- 
ues, correponding to the average distances between faults 
£ given in the caption. We checked that P{I)rs is unaf- 
fected for £ < TV/40. This shows that the extent of the 
four-point order has also an effect on the density proba- 
bility, much larger than in the Bernoullization procedure. 

Finally, let us mention that the problem of C- 
homometry, which is somehow the true homometry, is 
clearly related to the unicity of inverse problems, which 
is beyond the scope of this paper. It is important to 
note, however, that the use of ptychography [5], in which 
diffraction patterns are obtained by shifting illumination 
on the sample, can solve difficult problems of phase re- 
trieval [29|| , including C-homometry. 



DISCUSSION 

It might seem pointless to discuss the problem of ho- 
mometry while phase retrieval algorithms and ptychogra- 
phy can provide the full structural information, includ- 
ing high-order correlation functions. However, the full 
measurement of 3D speckle patterns is time consuming 
and in many situations it is not possible to get the data 
needed for such inversions. Methods of speckle analy- 
ses on quickly measured diagrams are thus needed to get 
novel information on the materials. 

As stressed in ref. [30] , real cases of true B-homometry 
are rare. On the contrary, because disorder is concerned, 
D-homometry is very frequent especially when systems 
are large. Moreover, it is related to high-order CF, which 
are hardly accessible to experiments [1, [Tq|, [31 ] 32]. How- 
ever, although D-homo metric systems have almost al- 
ways different speckle patterns, it is not yet clear whether 
this difference is quantifiable. Indeed, though the proba- 
blity density Prs (I) of our test bed sequence is strikingly 
different from its homo metric BS, the observed effects are 



very sensitive to disorder, and are almost invisible when 
four-point correlation lengths are too short. In this re- 
spect, we checked that the probablity density of the lat- 
tice shown in Fig. [I] does not present sizeable deviation 
from the normal law. Qualitatively, the "spikiness" of 
speckle patterns is quickly reinforced by the introduction 
of randomness, which makes deviation to the negative 
exponential curve delicate to observe. 

Consequently, though it is tempting to conclude from 
this study that high-order correlation functions are di- 
rectly observable through speckle statistics analysis, 
much theoretical work is still needed to find the rele- 
vant parameters controlling the statistics. Such an effort 
could be supported by more sophisticated analyses such 
as the use of second-order (or higher) probability func- 
tions [25]. The simple fact that the lattices of Fig. [T] 
can be reconstructed with minimum information shows 
that high-order CF are somehow hidden in the speckle 
repartition. 

An experimental difficulty lies in the presence of two- 
point correlations in all real systems (the classical SRO), 
which could obviously mask the speckle distribution anal- 
ysis. In this respect, we have checked that, at least for 
SRO lattices with no high-order correlations, dividing 
the speckles pattern intensity by its associated diffuse 
scattering one Idd (by smoothing, averaging or fitting) 
makes P(I/Idd) follow the normal decreasing exponen- 
tial law. This can help disentangling high-orders effects 
from two-point ones. 

Another issue might be the partial coherence of the 
beam, which reduces the speckles contrast and makes the 
previous analyses difficult. This could be overcome by 
the analysis of the speckles maximum intensities, which 
exhibit similar statistical properties (not shown here). 

In conclusion, we suggest the that speckle statistics 
analysis could be used as another tool to test for the 
presence of high-order correlations. Indeed much physics 
could be explored with the measurement of high-order 
correlations, still hidden to experiments (see ref. [3l| for 
examples). We hope this work will impulse theoretical 
and experimental studies on the role of high-order corre- 
lation functions on coherent diffraction diagram. 

We thank F. Berenguer, D. Gratias, D. Le Bolloc'h 
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